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G\ Abstract 

O ' 

C^i ' Here we give a brief overview of some investigations of the gravitational microlensing influence on broad absorption spectral 
lines of lensed QSOs. Especially, we consider the microlensing influence on the Fe Ka broad absorption lines using a model of an 
Q_! accretion disk covered by an absorption region. Gravitational microlensing is modeled by ray shooting method which enables us to 
D obtain realistic microlensing patterns. We obtain that microlensing can affect both emission and absorption component of line that 
depends on dimensions on emission and absorption line regions. Here we give detailed analysis of emission and absorption line 
0^ shape variations due to gravitational microlensing. 
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i 1. Introduction 

o ■ 

Around 10% of type I (broad emission line) Active Galac- 
&i tic Nuclei (AGN) sh ow gas outflowing from their centers 



C3 



dReichard et al. . 2003 ) with velocities ~ 10 3 ~ 4 km s _1 (obtained 
from blue-shifted absorption lines). Such outflo ws are proba- 
bly launched from the central engine of quasars (Elvis, 2000). 



The broad absorption lines (BALs) can be dete cted mainly in 



CN) * ne UV P art of spectra, but al so in Fe Ka line (iChartas et al 
t^j- l2004l2007trDone et allEo07l) . 

Gravitational microlensing has proven to be a powerful probe 
of the structure at the heart of quasars and therefore microlens- 
ing of Broad Line Absorption Regio ns (BALR) can be used 
<3\ to study th e BALR structure (see e . g. Lewis and Bella, 1 1998; 
O IChelouchel 120031: ILewis and Ibatal 120031: IChelouchd. Eool 



etc.). 
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The X-ray emission of AGN could be significantly absorbed 
by an outflowing wind, especially in the case of so-called Low 
Ionization Broad Absorption Line (LoBAL) quasars. Pres- 
ence of such X-ray absorbers is confirmed in gravitationally 
lensed LoBAL quas a r H1413+117 (Clove rleaf) at z = 2.56 



(Cha rtas et all 120071) . IChartas et alj ([2007) found that spectra 



of this quasar show emission features redward of the Fe Ka 
line rest-frame energy and broad absorption features blueward 
of this energy. The line is only significant in the brighter image 
A and a microlensing event could explain its energy-dependent 
magnification. A cluster of galaxies at redshift z — 1 .7 con- 
tributes to the lensing of this system, but the absorption from 
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the lens does not have any significant effect on spectral proper- 
ties ofH1413+117 while, at the same time, its intrinsic absorp- 
tion is much larger and it could sign i ficant ly affect the X-ray 



emission of this quasar (IChartas et al 
In several previo us pape rs 



20041). 



20031: IPopovicetal 



Jovanovic and Popovic, 



2006a; 
20071 



Jovanovic and Popovic , 



2009) 



dPopovic et all 12001 , 
IPopovic etall l2006b : 
Ijovanovic et all 12008 : 
we have shown that the ef- 
fect of strong gravitational lensing can have influence on the 
intensity and shape of Fe Ka line. Consequently, this effect 
can be used as a probe for the structure of the innermost part of 
lensed quasars. 

The focus of this work is to study how much the combined 
influence of intrinsic absorption and gravitational microlensing 
could change the Fe Ka spectral line profile, emitted from X- 
ray emitting region of LoBAL QSOs. The paper is organized as 
following: in §2 we describe how to model the X-ray emission 
of accretion disk of QSOs, in §3 and 4 the combined influence 
of intrinsic absorption and gravitational microlensing on this 
emission is given, in §5 we describe the adopted parameters of 
X-ray emitting and absorption regions, and also those which 
are used for generating the microlensing pattern for image A of 
H1413+1 17. In §6 we present the main results of our numerical 
simulations and finally in §7 we outline our conclusions. 



2. Modeling the X-ray emission of accretion disk 

The X-ray radiation of AGN is originating from the inner- 
most part of the accretion disk and the shape of the most promi- 
nent spectral line in this range, Fe Ka, strongly depends on 
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Figure 1: The modeled accretion disk (left) and the corresponding Fe Ka emission (right). 



emissivity law of the disk. It is usually accepted that sur- 
face emissivity of the disk varies with radius as a power law: 
e(r) = so ■ r ? , where eo is emissivity constant and q is emissivity 
index. We modeled the emission of an optically thick and geo- 
metrically thin accretion disk around a supermassive black hole 
(SMBH) using numerical simulations based on a ray-tracing 
method in Kerr metric, taking into account only photon tra- 
jectories reaching the observer's sky plane in the infinity. In 
this method one divides the image of the disk on the observer's 
sky into a number of small elements (pixels). For each pixel, 
the photon trajectory is traced backward from the observer by 
following the geodesies in Kerr space-time, until it crosses the 
plane of the disk. Then, flux density of radiation emitted by the 
disk at that point, as well as redshift factor of the photon are 
calculated. In that way, one can obtain color images of the ac- 
cretion disk (see Fig. [1] left) which a distant observer would see 
by a high resolution telescope. The simulated line profiles can 
be calculated taking into account the intensities and received 
photon energies of all pixels of the corresponding disk image 
(see F ig. Q] ri ght). For more details about this method see e.g. 
Jovanovic and Popovid (120091) and references therein. 



3. A simple model of Fe K« absorption region 

There are different models of the X-ray absorbing/obscuring 
regions, like absorbing medium c omprised of cold absorbing 
cloudlets by iFuerst and Wu (2004), b ut here we wil l focus on 
the model given by Ijovanovic and Popovid d2007l) . In this 
model, absorption region is considered to be composed of a 
number of individual spherical absorbing clouds with the same 
small radii, scattered in space so that projections of their cen- 
ters to the observer's sky plane (Xj, F,) have bivariate normal 
distribution A/2 (fi, 2). Here, fi = [ptx,fiy] T ar, d 



2 = 



" x 
po-xO-y 



pCTxCTY 

crl 



where fix and jiy are the means of Xi and T,, crx and cry are their 
standard deviations and p is the correlation between them. We 
use the following values of these parameters in our numerical 
simulations: p — 0, fix = X A , fly = Ya and crx — cry — Ra, 
where (Xa, Ya) is the center and Ra is the radius of projection 
of entire absorption region. 

The absorption coefficient A(X, Y) f or every spherical cloud 
in the absorption region is given by djovanovic and Popovid 
2007): 



A(X,Y) = (l-I A (X,Y))-e 



g(X,Y) E Q -E A 



(1) 



where absorption intensity coefficient I A (X, Y) describes the 
distribution of absorption over the whole region, Ea is the cen- 
tral energy of absorption and cr E is the width of absorption band 
(velocity dispersion). 

To compare (qualitatively) our simulated profiles (Fig. |2j 
right) with observ ations, let us recall the results obtained by 
Done et al.l (12007b . They found an evidence for so-called P 
Cygni profile of the Fe Ka line in narrow line Seyfert 1 galax- 
ies. According to these authors, complex X-ray spectra of these 
objects show strong "soft excess" below 2 keV and a sharp drop 
at ~ 7 keV which can be explained either by reflection or by ab- 
sorption from relativistic, partially ionized material close to the 
black hole. They showed that a sharp feature at ~ 7 keV re- 
sults from absorption/scattering/emission of iron Ka line in the 
wind. In the case of 1H 0707-495, this absorption feature can 
be satisfactorily fitted by the P Cygni profile, where emission 
component at ~ 5 keV is followed by a strong absorption com- 
ponent at ~ 7 keV jDone et aUl2007b . 

In our model, we are able to change density of absorbers, size 
of absorption region (or covering fraction of disk emission) and 
also position of absorption region in respect to disk center. As 
an example, in Fig. [2] (right) we show a deformation of the Fe 
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Figure 2: The modeled accretion disk covered by warm absorbing matter (left) and the corresponding Fe Ka absorption, presented by dashed line (right); the solid 
line on right panel presents the Fe Ka without absorption. 



Ka line profile due to absorption region that covers approaching 
side of the disk (Fig. [2] left). As one can see from Fig. [2] in 
that case a P Cygni profile of the Fe Ka is present, having an 
emission compo nent at 5 keV and absorption one at 7 keV, as it 
was reported bv lDone et al.l (|2007). 



the corresponding magnification map using nu merical simula- 



4. Gravitational microlensing by a random star field in the 
lens plane 

The most realistic approximation for gravitational microlens- 
ing amplification is so-called quadrupole microlens. This 
model is applied to obtain a spatial distribution of magnifica- 
tions in the source plane (where an accretion disk of AGN is 
located), prod uced by a random st ar field placed in the lens 
plane (see e.g. IWambsganss , 1992 ). Such spatial distribution 
of magnifications is called microlensing map, microlensing pat- 
tern or caustic network. If we consider a set of N compact ob- 
jects (e.g. stars) which are characterized by their positions x,- 
and their masses m,;, then norm alized lens equation is given by 
(IJovanovic and Popovid 12009): 



1 - k c + y 

1 - K r 



■ y 



(2) 



where x and y are normalized image and source positions, re- 
spectively. The sum describes light deflection by the stars and 
the last term is a quadrupole contribution from galaxy con- 
taining the stars, where k c is a smooth surface mass density 
and y is an external shear. The total surface mass density or 
convergence can be written as k — k, + k c , where k„ repre- 
sents the contribution from the compact microlenses. The cor- 
responding microlensing map is then defined by two param- 
eters: the convergence - k, and the shear due to the external 
mass - y. For some specific microlensing event one can model 



tions based on ray-shoo ti ng techniques (see e.g. Popovic et ah 
2006b; IJovanovic et all 120081: IJovanovic and Popovid 120091 



and references therein), in which the rays are shot from the ob- 
server to the source, through a randomly generated star field in 
the lens plane. These light rays are then collected in pixels in 
the source plane, and the number of rays in one pixel is pro- 
portional to the magnification due to microlensing at this point 
in the source plane. As an example, the magnification map of 
Cloverleaf quasar image A is shown in Fig. [3] 




Figure 3: Magnification map of Cloverleaf quasar H1413+117 image A. A 
white solid line represents the analyzed path of an accretion disk center. 
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5. Combined influence of intrinsic absorption and gravita- 
tional microlensing on the Fe Ka emission of QSOs 

To obtain modeled line profiles, it is necessary to define a 
number of parameters which describe the line emitting region 
in the disk, such as constraints for its size, the disk inclination 
angle and the angular momentum of the central BH. For that 
purpose we usually use the results from stu dies of the observed 
Fe K a line profiles in AGN, such as e.g. iNandra et alJ (11997 , 
20071) . The inner radius /?;„ of the disk cannot be smaller than 
the radius of the marginally stable orbit R ms , that corresponds 
to R m <, = 6 R g (gravitational radius R g = GM/c 2 , where G is the 
gravitational constant, M is the mass of central BH, and c is the 
speed of light) in Schwarzschild metric and to R ms = 1 .23 R g 
in the case of Kerr metric with angular momentum parameter 
a = 0.998. To select the outer radius R out of the disk, we take 
into account some recent investigations of the Fe Ka line pro- 
file showing that it should be emitted from the innermost part 
of th e disk which outer ra dius is within several tens of R g (see 
e.g. IPopovic et all l2006al and references therein). Therefore, 
in our numerical simulations we adopted the following disk pa- 
rameters: inclination i - 35°, inner and outer radii 7?,„ = R ms 
and R ou , - 20 R g , emissivity index q = -2.5 and angular mo- 
mentum of central BH a — (Schwarzschild metric). 

In order to study the influence of intrinsic absorption of 
Cloverleaf quasar H1413+117 on its Fe Ka line profile, we 
modeled an X-ray absorption region using the following param- 
eters: ya - 0.2 R g « Ra -1 R g (and thus I A (X, Y) can be as- 
sumed as constant), central energy of absorption Ea = Eq = 6.4 
keV and velocity dispersion cr E = 0.5 keV (that corresponds to 
random velocity of clouds of ~ 20000 km/s). We assumed that 
the absorption region is located at the following positions: (i) 



-15 R g , Y A = QRg, (ii)X A 



R g , Y A 



R„ and (iii) 



X A = 15 R g , Y A = R g and that it consists of 10000 individual 
absorbing clouds. 

Microlensing magnification pattern for the image A of 
H1413+117 (see Fig. EJ with 16 Einstein Ring Radii (R E ) on 
a side (where Re ~ 1431 R g ) is calculate d using the followin g 
parameters: red shift of the lens: Zd — 1 .7 (Char tas et all 2004), 
convergence k — 0.52, shear y = 0.14 dKavseret al.Lll990h and 
the masses of individual deflectors (stars) mj = 1 M . Mi- 
crolensing time scales, which are necessary for obtaining the 
simulated light curves, are estimated from correspondin g dis- 



tance scales according to the formula (13) of iJovanovic et al 
(2008), in which /? SOU ice is replaced by the distance from the 



center of accretion disk, assuming a flat cosmological model 
with Q. = 0.3 and H = 75 km s^Mpc" 1 . 



6. Results 

The results of numerical simulations of a relativistic accre- 
tion disk partially covered by an X-ray absorption region where 
the Fe Ka line radiation is magnified due to gravitation mi- 
crolensing are presented in Fig. [4] (left) for three different po- 
sitions of absorption region with respect to the accretion disk 
center. In all three panels of Fig. |4] this system consisting of 



accretion disk and absorption region is located at the same po- 
sition on the microlensing map ofH1413+117 image A. The 
corresponding comparisons between the unabsorbed, absorbed, 
magnified and combined absorbed-magnified Fe Ka spectral 
line profiles are given in Fig. (right). As one can see from 
Fig. HI when X-ray radiation from approaching side of the disk 
is significantly absorbed there is a very strong absorption com- 
ponent of the iron line and in such case the emission Fe Ka 
component looks redshifted at ~ 5 keV and is followed by a 
strong absorption component at ~ 7 keV, indicating the P Cygni 
profile of the iron line (see also Fig. 

Gravitational microlensing affects both components, but the 
magnification of absorption component is much stronger (see 
dotted line on the top panel of Fig. @J. 

When the absorption region is located over the central part of 
the disk, it also causes the occurrence of P Cygni profile of the 
iron line, but in this case the intensity of absorption component 
is much smaller, even when we take into account magnification 
due to microlensing (see the middle panel of Fig. 0). 

Finally, in the third case, when the absorption region is lo- 
cated over the receding side of the disk, the absorption com- 
ponent of the iron line is missing and the line profile is just 
slightly changed due to intrinsic absorption. In the combination 
with gravitational microlensing influence, it causes significant 
deformations of the Fe Ka line profile, but also without any in- 
dication for the absorption component of the iron line (see the 
bottom panel of Fig. @J. 

The corresponding simulated light curves, produced when 
the system consisting of the accretion disk absorbed emission 
crosses over a magnification pattern along the horizontal path 
in Fig. [3] are given in Fig. [5] 

The three panels in Fig. [5]show the simulated light curves for 
three different positions of absorption region across the accre- 
tion disk. We compared simulated light curves in the case when 
the pure Fe Ka emission is magnified (solid line) and when ab- 
sorbed emission of the Fe Ka line is magnified (dashed line). 
As it can be seen from Fig. [5] there is a global correlation be- 
tween these two light curves. The greatest difference between 
the two light curves is in the case when absorption region cov- 
ers the central part of the disk and the smallest one is when it 
covers the receding side of the disk. Such behavior of the light 
curves is a consequence of different absorption rate of the iron 
line for three different positions of absorption region. 

From Fig. |5]we can also estimate the frequency of so-called 
high magnification events (HMEs - asymmetric peaks in the 
light curves which depend not only on microlens parameters, 
but also on parameters of the X-ray emitting region), i.e. the 
number of such events per unit time. According to the obtained 
results, in the case of image A of H1413+117 one can expect 
one such event approximately every 15 years, along the spec- 
ified path. But, there are also areas in the microlensing map 
of H1413 + 117A with a significantly higher number of caus- 
tics (see e.g. top part of Fig. [3]), and therefore a source tran- 
sit over these regions would produce HMEs with much higher 
frequency. In principle, it can be expected that the major- 
ity of HMEs should be detected in X-ray light curves, less of 
them in UV and the smallest number in optical light curves 
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Figure 4: Left: Relativistic accretion disk partially covered by an X-ray absorption region (randomly scattered gray dots). Magnification due to gravitation 
microlensing is also presented by contour lines. Positions of absorption region are: Xa = -15 R g , Ya = R g (top), Xa = R g , Ya = R g (middle) and Xa = 15 R g , 
Ya = R g (bottom). Right: Comparison between the following corresponding Fe Ka spectral line profiles: unabsorbed (solid line), absorbed (dashed line), 
magnified unabsorbed (dash-dot line) and absorbed and magnified profile (dotted line). 
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Figure 5: Comparison between simulated light curves of unabsorbed (solid line) and absorbed Fe Ka line flux (dotted line). Simulations are performed for the 
accretion disk crossing along the horizontal path in the magnification map of image A of Cloverleaf quasar H1413+1 17 (see Fig. |3j- The positions of absorption 
region with respect to the accretion disk center are taken as: Xa = — 15 R g , Ya = R g (top), X& = R g , = R g (middle) and Xa = 15 R g , Ya = R g (bottom). 
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(IJovanovic et all 120081) . IJovanovic et a l. (2008) found that the 



typical microlensing time scales for the X-ray band are on order 
of several months, while for the UV/optical bands they are on 
order of several years. 

7. Conclusions 

In order to study how much the combined influence of ab- 
sorption and gravitational microlensing could change the Fe 
Ka spectral line profile, emitted from the X-ray emitting region 
of LoBAL QSOs we developed a simple model of the X-ray 
absorbing/obscuring region that could explain the observed P 



Cygni profile of the Fe Ka line recently observed by lDone et al 



(2007). We found that the model can reproduce P Cygni profile 
only if at least a part of approaching side of the disk is covered 
by absorbing material (see Fig. [2}. In particular case of 1H 
0707-495, the observed sharp drop at ~ 7 keV and emission at 
5 keV can be explained by absorption from relativistic, partially 
ionized material close to the black hole that meanly cover the 
approaching side of the relativistic accretion disk. 

It is interesting that the width and depth of the absorption 
component strongly depends on the projection of this region on 
the accretion disk. 

On the other side, this model (absorbed emission of the 
accretion disk) has been used for simulation of gravitational 
microlensing of such system by a random star field in the 
lens plane. We estimated the magnification of the line flux 
and calculated the simulated light curves produced when such 
source crosses over the H1413+1 17A microlensing magnifica- 
tion map. Our results show that combined influence of absorp- 
tion and gravitational microlensing could explain the observed 
variations of the Fe Ka line profile in the case of H1413+1 17 
(emission features redward of the line rest-frame energy and 
broad absorption fe atures blueward of this energy) reported by 
IChartas etal](l2004 . 

Note here, that we made a net of models for variations of the 
Fe Ka absorption of H1413+117A due to microlensing which 
can be used for comparison with observed variations. We hope 
that these results (and model) will be used in order to probe 
both, absorption region and accretion disk characteristics in 
lensed QSOs with Fe Ka absorption component. 
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